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An estimate of the NNLO supersymmetric QCD eflects for Higgs production at iiadron colliders 
is given. Assuming an effective gluon-Higgs interaction, these corrections enter only in terms of 
process-independent, factorizable terms. We argue that the current knowledge of these terms up 
to NLO is sufficient to derive the NNLO hadronic cross section within the limitations of the stan- 
dard theoretical uncertainties arising mainly from renormalization and factorization scale variations. 
The SUSY contributions are small with respect to the QCD effects, which means that the NNLO 
corrections to Higgs production are very similar in the SM and the MSSM. 
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I. INTRODUCTION 

Gluon fusion is the dominant production mechanism for Higgs bosons at the Large Hadron Colhder (LHC) 
(for reviews, see Refs. 1, 2]). A feature of the gluon fusion process is that it is loop- mediated already at 
leading order. This makes it particularly sensitive to non-standard particles and couplings as they are 
predicted by extended theories. A very popular extension of the Standard Model (SM) is the Minimal 
Supersymmetric Standard Model (MSSM) (for a review, see Ref. Q), with its five physical Higgs bosons. 

We will focus on a scenario where the ratio of the vacuum expectation values of the two Higgs doublets 
is not too large, tan/3 ^ Mt/Mi,, (Mt=top mass, Mb = bottom mass), so that the bottom is much 
smaller than the top Yukawa coupling. In this case, the dominant effects on the gluon-Higgs coupling 
in the MSSM arise from the top quark t and its scalar supersymmetric (SUSY) partner, the top squark 
i. SUSY-QCD corrections are induced by virtualgluons g and their fermionic SUSY partners, the gluinos 
g. These effects have recently (see also Ref. U) been evaluated at next-to-leading order (NLO) in the 
limit where ^ {Mt, Mj, Mg}, where (f> denotes either of the two CP-even Higgs bosons, h or H. This 
limit is expected to work extremely well, if the leading order (LO) dependence on Mt, Mj:, Alg is taken 
into account exactly. This can be inferred from the NLO behavior in the SM 0, 0, B 13 ■ 
In the effective Lagrangian approach, the evaluation of the hadronic Higgs cross section factorizes into the 
calculation of the effective gluon-Higgs coupling, times the calculation of the actual process 4>+X as 
mediated by the effective gluon-Higgs operator. For a full next-to-next-to-leading order (NNLO) result in 
this approach, both factors need to be evaluated up to NNLO. However, in the SM, the NNLO contribution 
of the effective coupling leads to a numerically negligible contribution, and we will argue that this is true 
also in the MSSM. The NNLO Higgs production cross section can therefore be evaluated from the NLO 
expression of the effective coupling, as taken from Ref. and the NNLO results for the process diagrams, 
which are identical to the SM case 0, ^j, . 
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II. THE APPROXIMATION 



A. Definition and Standard Model case 



We use the effective Lagrangian approach where the top quark and all supersymmetric particles are 
considered heavy with respect to the Higgs boson, see Ref. 4]. In this case, the hadronic cross section 
ffhk = (j{hk + X) for Higgs production can be written as 

cThk{z) ^ (ToC^'^hkiz) , (1) 



j ■ Jz Jz/xi 
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where i, j denote any partons inside the hadrons ft., fc, and ipi/h{x) are the parton densities; is the 
Higgs boson mass, and s is the hadronic center-of-mass (cm.) energy. 

The coefficient function C, defined below, contains the remnant dependence of the gluon-Higgs coupling 
on the heavy masses, and co is defined such that the leading order dependence on these masses of ahk{z) is 
exact. Its exact form is irrelevant for our argument and shall not be given here, owing to space limitations 
(see, e.g. Ref. 0). 

The partonic expression can be expanded in terms of Ug: 



E., {x) = {x) + ^ ^\\> {x) + [^) T,Y> (x) + 0{at) 



.(2), 



(3) 



where x = M^/s, and s is the partonic cm. energy. Here and in what follows, ag denotes the MS- 
renormalized strong coupling constant for five active quark flavors, evaluated at the renormalization scale 
MR- 

For the coefficient function, we write 



(4) 



In the SM, the Hi are known for i = 1, . . . , 3 00 (K3 contributes only at N^LO). In the MSSM, ki has 
been evaluated only recently jjj. 

We now define 



4^w=E r^xi r dx, 

i i J z J z/xi 



V%/h{xi) ^j/k{x2) S 
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(5) 



For the s|;^.^(z), ri = 0,1,2, we assume that the parton densities fi/h are evaluated at NNLO.^ Thus, the 
NNLO expression for the hadronic cross section can be written as 



„NNLO 

CT =0-0 



S(0) + ^ fs(i)+2«il]( 

TT V 

(e(2)+2«iI](i) + («2 + 2«2)S 



(6) 



^ We use the approximate NNLO parton densities of Ref. Il6l . 
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FIG. 1: Individual NNLO contributions to the total hadronic Higgs production cross section. The 
notation is defined in Eqs. and ||HJ . The renormahzation and factorization scale fin and /ip are 
identified with the Higgs mass M^. 
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where the indices h,k E have been dropped for simphcity. The basis of our estimate of the NNLO 

terms in SUSY will be that the numerical effect of the term proportional to K2 in Eq. lO is negligible 
compared to the theoretical uncertainty of the NNLO prediction. 

For that purpose, let us look at the relative magnitude of the S'"-' in the case of Higgs production at 
the LHC, Fig.n(a). We see that E^"^ is more than two orders of magnitude smaller than which 
suggests that the effects from H2 can be neglected, if K2 is not too large. In order to get a feehng for the 
magnitude of K2, let us look at the SM case. There we have 

^(o),SM^_l^ ^SM^ 11 ^2.75, 

^ ^ (''') 

SM^ 2777 19 (Jl^\i\ "ij' 6.153 + 2.854 Z.t, 

2 288 16 V 96 3 ^ ) ' 

with l^t = ln-(/i^/M(^), where is the renormahzation scale and Mt is the on-shell top quark mass. 

Using these numbers, we arrive at Fig.^(6). It shows the relative size of the four terms that contribute 
to the cross section in Eq. © at order a*: 

Xi = S(2), X2 = 2«:iSW, X3 = «?SW, X4^2k2^^°'>. (8) 

As expected from the numerical value of k^^, Eq. ||7J), X4 is indeed negligible with respect to Xi : it is 
down by a factor of 30. But another remarkable observation is that the term proportional to E^^-*, i.e. 
X2, amounts to around 30% of the full contribution. For comparison: the {2 ki E^^^) term in Eq. ||SJ) 
amounts to only 15% of the complete contribution. 

To summarize: In the SM, the term K2 to the coefhcient function of Eq. gives a negligible con- 
tribution to the NNLO cross section. In fact, we checked that the difference between the true^ and the 
approximate NNLO cross section (i.e. with K2 ~ 0) is less than 1% at the LHC. This is much smaller 
than the theoretical uncertainty of around 15%, as estimated by the variation of the factorization and 
the renormahzation scale at NNLO. On the other hand, the knowledge of ki is, relatively speaking, nu- 
merically more important for the NNLO than for the NLO contribution to the cross section, for which it 
was originally evaluated Q . 



^ Within the eflfective theory approach. 
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FIG. 2: K-factor in the Standard Model (dashed blue) and the MSSM (solid red) for the indicated set 
of parameters. The narrow (red) band in the NNLO MSSM curve corresponds to varying K2 between 
zero and 2k2^. The renormalization and factorization scales {^n, Mf) have been identified with the 
Higgs mass in these curves. The diagonally shaded band (green) corresponds to the variation of /ir 
between 2M^ and in the NNLO result with K2 = nf^ (fiF = M^). 



B. Minimal Supersymmetric Standard Model 



In the MSSM, we can parameterize the NLO corrections to the effective Lagrangian as 

^SUSY^^SM^^^^^ n^^^^^ (9) 

In addition, the tree-level normahzation of Eq.jBjl. C^^\ changes, of course, but this is irrelevant for our 
discussion. Sni was recently computed in Ref. jj]"^ and was shown to be negative, with 

l'5^i|<l (10) 

for relevant values of the SUSY parameters (recall that we restrict ourselves to tan/3 <^ Mt/Mt,). It is 
thus reasonable to assume that also the value of K2 in SUSY-QCD will be of the same order of magnitude 
as in the SM (or smaller) . Combining this assumption with the discussion of Fig.^ (see above) leads us 
to the conclusion that the NNLO cross section for hadronic Higgs production in supersymmetry should 
be approximated well by setting K2 ~ k^^^. 



III. RESULTS 



As in Ref. , we will neglect squark mixing and set the bottom Yukawa coupling to zero for simplicity. 
More detailed phenomenological studies have to be deferred to a forthcoming publication. 

Fig.Elshows the NLO and the NNLO K-factor, A'x = o-^/cr^° (X = NLO, NNLO) in the SM case (dashed), 
and in the MSSM, for Mf = Mt = 175 GeV, and Mg = 500 GeV; cr^o, (7^^°, and (jNNLO ^^^.g evaluated 
with LO, NLO, and NNLO parton densities and as evolution. 

The NNLO result in the MSSM is given by the narrow (red) band, arising from the variation of K2 between 
zero and 2 k™ (see Eq. (0). This should serve as an estimate of the theoretical uncertainty induced by 
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the approximation introduced in Sect.^ Within our approximations, the K-factor in Fig.|21is vahd for 
both CP-even Higgs bosons of the MSSM, since the Yukawa couphng cancels in the ratio of the LO to the 
higher order resuhs. 

As at NLO, the SUSY effects are small with respect to the QCD effects at NNLO, so that the total K- 
factor in the MSSM is very similar to its SM value. The theoretical uncertainties due to variation of 
the renormalization and factorization scales are almost identical to the SM case, since the only source of 
additional scale dependence in the MSSM arises from K2- For K2 = the /iR dependence of the NNLO 
prediction is indicated in Fig.|21as the diagonally shaded band; the fj.p dependence is much smaller. 



IV. CONCLUSIONS 



We have discussed an approximation for the NNLO contributions to supersymmetric Higgs production in 
gluon fusion. It was argued that the corresponding NNLO estimate should be accurate up to a few per 
cent, which is much smaller than the theoretical uncertainty induced by the residual renormalization and 
factorization scale dependence, as well as the anticipated experimental accuracy. For a given set of SUSY 
parameters, and within the restrictions on these parameters as discussed in the main text, the production 
cross section for a Higgs boson in the MSSM is thus known to a precision similar to that in the SM. More 
detailed studies of the MSSM parameter space are clearly desirable and will be presented elsewhere. 
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